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The successful operation of X-ray free-electron lasers (FELs), like the Linac Coherent Light Source 
or the Free-Electron Laser in Hamburg (FLASH), makes unprecedented research on matter at atomic 
length and ultrafast time scales possible. However, in order to take advantage of these unique light 
sources and to meet the strict requirements of many experiments in photon science, FEL photon 
pulse durations need to be known and tunable. This can be achieved by controlling the FEL driving 
electron beams, and high-resolution longitudinal electron beam diagnostics can be utilized to provide 
constraints on the expected FEL photon pulse durations. In this paper, we present comparative 
measurements of soft X-ray pulse durations and electron bunch lengths at FLASH. The soft X-ray 
pulse durations were measured by FEL radiation pulse energy statistics and compared to electron 
bunch lengths determined by frequency-domain spectroscopy of coherent transition radiation in 
the terahertz range and time-domain longitudinal phase space measurements. The experimental 
results, theoretical considerations, and simulations show that high-resolution longitudinal electron 
beam diagnostics provide reasonable constraints on the expected FEL photon pulse durations. In 
addition, we demonstrated the generation of soft X-ray pulses with durations below 50 fs (FWHM) 
after the implementation of the new uniform electron bunch compression scheme used at FLASH. 

PACS numbers: 29.27.-a, 41.60.Cr, 41.50.+h 



I. INTRODUCTION 

Since the first operation of a laser in 1960 [T], a tremen- 
dous number of new experimental techniques became 
possible with continuously changing and growing require- 
ments on the laser systems, e.g., higher spectral bright- 
ness, shorter wavelengths, or photon pulse durations with 
simultaneous tunability. High-gain free-electron lasers 
(FELs) meet many of theses requirements, and the suc- 
cessful operation of X-ray FELs like the Free-Electron 
Laser in Hamburg (FLASH) J] , the Linac Coherent Light 
Source (LCLS) j3], or the SPring-8 Angstrom Compact 
Free Electron Laser [3] make unprecedented research on 
matter at atomic length scales possible [5], The demon- 
stration of FEL photon pulse durations in the femtosec- 
ond range (e.g., Refs. [HIE]) has further extended the ca- 
pabilities to research of dynamical processes at ultrafast 
time scales (see, e.g., Ref. [5]) with ongoing demands on 
the generation and control of ultrashort photon pulses. 

In recent years, several methods to control FEL pho- 
ton pulse durations have been proposed by manipulat- 
ing and controlling the FEL driving electron bunches. 
The low-charge operation at LCLS [9J demonstrated elec- 
tron bunch lengths below 10 fs [TO], and the same strat- 
egy of low-charge operation is planned for the European 
XFEL and FLASH 11J. Other methods, with additional 
prospects of generating photon pulses in the attosecond 
range, make use of electron bunch manipulation with 
conventional quantum lasers (e.g., Refs. P^HTfJ] ) or by 
selectively spoiling the transverse emittance of the elec- 
tron beam [17]. However, reliable operation of theses 
methods requires capabilities to diagnose the FEL pho- 



ton pulse shapes and durations with high accuracy, which 
is a tremendous challenge and an active field of research. 

First single-shot measurements of FEL photon pulse 
durations with femtosecond accuracy have recently been 
demonstrated by terahertz-field driven streaking exper- 
iments in the time-domain [T5] ITS"] , and the statistical 
and spectral properties of FEL radiation emitted in the 
exponential gain regime allow to measure the mean pho- 
ton pulse duration (e.g., Refs. P0"H25] ). Other proposed 
methods make use of the FEL induced slice energy loss 
in order to measure "replicas" of the FEL photon pulses 
(see Refs. [Ml 123). As is presented in this paper, high- 
resolution longitudinal electron beam diagnostics in stan- 
dard configurations, i.e., as commonly used at present 
X-ray FELs, can provide reasonable and complementary 
constraints on the expected FEL photon pulse durations. 

In this paper, we compare measurements of electron 
bunch lengths with corresponding soft X-ray pulse dura- 
tions at FLASH, and discuss their relationship by theo- 
retical considerations and simulations. The experimental 
setup and applied methods are described in Sec. In] and 



in Sec. Ill we discuss the experimental results for FEL op- 



eration with different bunch charges, which corresponds 
to different electron bunch lengths and soft X-ray pulse 
durations. The final results are compared in Sec. [TV 
the summary and conclusions arc presented in Sec 



II. EXPERIMENTAL SETUP AND METHODS 

The measurements presented in this paper were car- 
ried out at FLASH, which is a self- amplified spontaneous 
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FIG. 1: Layout of the Free-Electron Laser in Hamburg (FLASH) with its superconducting (SC) accelerating structures (ACC), 
the two magnetic bunch compressor (BC) chicanes, and the new third-harmonic rf linearizer system (L3). The positions of the 
diagnostics used for photon pulse duration and electron bunch length measurements are indicated by green dots and arrows. 



emission (SASE) FEL for extreme- ultraviolet and 
soft X-ray radiation, driven by a superconducting radio- 
frequency (rf) linear accelerator [12] . The schematic lay- 
out of FLASH is depicted in Fig. [T] showing the injec- 
tor, which is based on a laser-driven normal conducting 
rf gun, the superconducting accelerating structures, two 
magnetic bunch compressor chicanes, and the undulator 
magnet system. The positions of the diagnostics used for 
FEL photon pulse duration and electron bunch length 
measurements are indicated by green dots and arrows. 

In 2010, FLASH was upgraded [37] to a maximum 
beam energy of 1.25 GeV, allowing for the generation of 
soft X-rays below 4.4 nm (water window) with its funda- 
mental. One of the major upgrades was the installation of 
the third-harmonic rf system, which is dedicated for the 
linearization of the longitudinal phase space upstream of 
the first bunch compressor [28l|29]. Before this upgrade, 
FLASH was operated with non-uniformly compressed 
electron bunches with a short leading spike [3TJ1 I5T| . 
The corresponding FEL photon pulses had short dura- 
tions but low pulse energies due to the fact that only a 
small fraction of the bunch (the short leading spike) con- 
tributed to lasing [2J[B]. In Sec. Ill it is shown that the 
third-harmonic rf linearizer system permits more uniform 
bunch compression with higher FEL photon pulse ener- 
gies and simultaneous tunability of the pulse durations. 

In the following, we describe the experimental setup 
and applied methods to measure the FEL photon pulse 
durations and electron bunch lengths. The photon pulse 
durations were estimated by the statistics of FEL radia- 
tion pulse energies measured with a micro-channel plate 
based detector. The electron bunch lengths were mea- 
sured in the time-domain by using a transverse deflecting 
rf structure in combination with a magnetic energy spec- 
trometer and in the frequency-domain by spectroscopy 
of coherent transition radiation in the terahertz range. 



A. Time-domain longitudinal phase space 
diagnostics for electron beams 

The time-domain electron bunch length measurements 
were done by using a transverse deflecting rf structure 
(TDS) [32^ in combination with a magnetic energy spec- 
trometer. Transverse deflecting rf structures are widely 
used for electron bunch length and longitudinal pro- 
file measurements at present FELs, and provide high- 



resolution single-shot diagnostics (see, e.g., Ref. [3~3"]). 
In combination with standard techniques for transverse 
emittance measurements, the operation of TDSs makes 
time-resolved emittance, i.e., slice emittance, measure- 
ments possible. The complementary use of an energy 
spectrometer allows direct longitudinal phase space and 
slice energy spread measurements (e.g., Refs. [3T] |3"H - 
136] L Recently, a TDS in combination with a magnetic 
energy spectrometer downstream of FEL undulators have 
been proposed to measure the FEL-induced slice energy 
spread for temporal X-ray pulse characterization |25j . 

At FLASH, a LOLA-type [32] TDS was successfully 
operated upstream of the energy collimator [31] before 
it was moved close to the undulators during the FLASH 
upgrade in 2010 [57] . The TDS has been integrated in 
a dedicated setup for measurements of the longitudinal 
phase space (35] [36] . As depicted in Fig. [I] the TDS can 
either be operated in combination with the dispersive 
energy spectrometer or by using off-axis screens in the 
non-dispersive main beamline during FEL operation. In 
both cases, the screen stations are equipped with differ- 
ent imaging screens and a camera system with motorized 
optics. A fast kicker magnet (not shown in Fig. [IJ can 
operate the off-axis screens for TDS measurements and 
frequency-domain electron bunch length diagnostics by 
terahertz spectroscopy (see Sec. II B I in parallel. Tech- 



nical details and performance measurements on the new 
measurement setup can be found in Refs. [33 ES], and 
technical information about the TDS at FLASH and de- 
tailed descriptions of time-domain electron bunch diag- 
nostics using a TDS can be found in Refs. [SUET]- Here, 
we describe only the basic principles of longitudinal phase 
space diagnostics required throughout this paper. 

The vertical betatron motion of an electron passing a 
TDS around rf phase zero-crossing is given by [311 137] 



U±(s) = 2/o (s) + C y (s,s )c x z ± S y (s, s )c 



(1) 



with the linear correlation C y (s,so) and shear function 
S y (s,s ) = R 34 K y = v //3 y ( S )/3. y ( So )sin(A0 y ) e ^ , (2) 



where R34 — f3 y (s) (3 y (s )sm(A(f>y) is the angular-to- 
spatial element of the vertical beam transfer matrix from 
the TDS at so to any position s, (3 y is the vertical beta 
function, Acj) y is the vertical phase advance between sq 
and s, and j/o describes the intrinsic offset. The expres- 
sion K y = eujVy / E is the vertical kick strength with the 
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peak deflection voltage V y in the TDS, c is the speed 
of light in vacuum, e is the elementary charge, E is the 
electron energy, z is the longitudinal position of the elec- 
tron relative to the zero-crossing rf phase, and u)/(2tt) is 
the operating rf frequency. The longitudinal-to-vertical 
correlation C y in Eq. ([I]) is independent of the TDS op- 
eration and can exist intrinsically due to time-dependent 
kicks generated from collective effects such as coherent 
synchrotron radiation or wakefields. This linear correla- 
tion may lead to systematical errors in electron bunch 
length measurements, which, however, can be removed 
by performing measurements at two zero-crossing TDS 
rf phases shifted by 180 deg, i.e., with ±S y (cf. Eq. dTl). 

The expression in Eq. (JlJ shows a linear mapping 
from the longitudinal to the vertical coordinate and al- 
lows longitudinal electron bunch profile measurements 
by means of transverse beam diagnostics using imaging 
screens or wire-scanners. The vertical shear function S y 
determines the slope of this mapping and can be cali- 
brated by measuring the vertical centroid offset of the 
bunch as a function of the TDS rf phase. The electron 
bunch current is given by the calibrated and normalized 
longitudinal bunch profile multiplied by the measured 
electron bunch charge. In the following, the longitudi- 
nal coordinate z is expressed by a time coordinate via 
t = —z/c and the bunch length by the corresponding 
bunch duration. The latter is expressed by the root mean 
square (r.m.s.) value E t>e = \C y ± S' v | _1 (c^ ± - cr^) 1 / 2 , 
where a y ± is the measured vertical r.m.s. beam size dur- 
ing TDS operation corresponding to ±S y , respectively. 
The r.m.s. bunch duration S t describes a quadratic 
equation for cr 2 ± , and by performing measurements of 
the vertical r.m.s. beam size at ±S y and S y = (TDS 
switched off), the longitudinal-to- vertical correlation C y 
and the intrinsic vertical r.m.s. beam size a y o can be 
determined. These parameters limit the r.m.s. time res- 
olution to Ht, e = v y o/\C y ± S y \ (cf. Refs. [SUEZ]). 

By combining the operation of a TDS with an en- 
ergy spectrometer and using imaging screens to get two- 
dimensional transverse beam profiles, longitudinal phase 
space (energy versus time) measurements with single- 
shot capability can be achieved. The simplest magnetic 
energy spectrometer consists of a dispersive beamlinc 
downstream of a dipole magnet. The corresponding hor- 
izontal betatron motion, which is perpendicular to the 
shearing plane of the TDS, is given by (e.g., Refs. [3TII3?] ) 

x(s) =x (s) + D x (s, s )S (3) 

with the horizontal dispersion D x and the relative mo- 
mentum deviation 5 — Ap/p. The magnetic energy spec- 
trometer used at FLASH has a nominal horizontal dis- 
persion of about 750 mm at the position of the screen 
station [33]. For relativistic electron beam energies with 
a Lorentz factor of 7 3> 1, which is the case throughout 
this paper, the electron beam energy is given by E pc, 
and S can be described as the relative energy deviation. 
Then the expression in Eq. ^ represents a linear map- 
ping between the relative energy deviation and the hori- 



zontal coordinate, where D x determines the slope of the 
mapping. The dispersion D x can be calibrated by mea- 
suring the horizontal centroid offset of the bunch as a 
function of the energy deviation. The r.m.s. resolution of 
the relative energy deviation is defined as TZs.c — &xQ /D x 
(cf. time resolution lZt,c), but here, the intrinsic beam 
size, and correspondingly the energy resolution cannot 
simply be measured by switching the dipole magnet off. 
A possibility to estimate the relative energy resolution 
is given by the deviations in the measurement of a well- 
known energy spread (see Refs. [3TJ|3J}] for more details). 

B. Frequency-domain terahertz spectroscopy of 
coherent transition radiation 

The frequency-domain longitudinal electron bunch 
profile measurements were carried out by spectroscopy 
of coherent transition radiation using a multi-channel 
terahertz (THz) spectrometer. The schematic and ba- 
sic layout of the THz spectrometer is shown in Fig. [2] 
The design is based on five stages of blazed reflection 
gratings in combination with focusing ring mirrors. Fast 
readout of the pyroelectric detector arrays allows single- 
shot diagnostics. The technical details about the focusing 
ring mirror design and the spectrometer principle can be 
found in Ref . |38j , and technical details and performance 
measurements on the particular THz spectrometer used 
throughout the measurements presented in this paper can 
be found in Ref. [33]. At FLASH, two identical THz spec- 
trometers exist, one upstream of the energy collimator 
(see Fig. [T]) , which was used for the spectroscopic mea- 
surements presented in this paper, and one close to the 
undulators, which is currently being commissioned. 

The basic principle of spectroscopy for longitudinal 
electron bunch profile diagnostics relies on the spec- 
tral intensity of transition radiation (diffraction or syn- 
chrotron radiation are also suitable) emitted by an elec- 
tron bunch with TV electrons which is given by [39H42] 

% = d ^(N+(N-l)N\F(\)f) , (4) 

where dZ± /dX is the spectral intensity of a single electron, 
and F(X) is the longitudinal form factor of the bunch as 
function of the wavelength A, which can be expressed 
by the Fourier transform of the normalized longitudinal 
charge density pn{z), i.e., the electron bunch current, as 

/oo 
p N (z)e- i2wz / x dz. (5) 
-OO 

In Eq. @ we have made use of the fact that radiation 
from relativistic electrons is confined to small angles, 
hence we have neglected transverse beam size effects [33] • 
Inverse Fourier transform permits the reconstruction of 
the longitudinal electron bunch profile, i.e. the bunch 
current I(z) = Qpn(z) with the bunch charge Q, when 
the complex form factor F(X) — |.F(A)|e**( ' is known. 



4 




as demonstrated in Sec. IIIB and uncertainties due to 
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(b) One stage with grating, ring mirror, and detector. 

FIG. 2: Principle of the multi-channel THz and infrared spec- 
trometer based on five stages of blazed reflection gratings in 
combination with focusing ring mirrors. The schematic layout 
is depicted in (a), and (b) shows one stage exemplarily (see 
Refs. |38l I39| for technical details). The fast readout of the 
pyroelectric detector arrays allows single-shot spectroscopy. 



Only the modulus of the longitudinal form factor 
|F(A)| can be determined experimentally (see Eq. Q), 
and the phase information 3(A) remains unknown. How- 
ever, a phase retrieval can be achieved by applying the 
Kramcrs-Kronig relations (e.g., Refs. [30H12])) which con- 
nect the real and imaginary parts of a complex func- 
tion. Taking the logarithm of the complex form factor 
ln_F = In \F\ +i$ and applying a Hilbert transform yield 
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where 5>b is the Blaschke phase, V denotes the Cauchy 
principal value, and u) — 2nc/X. The Blaschke phase can- 
not be determined from the modulus of the form factor 
and is omitted ($b(w) = 0) in the following. A profile 
reconstruction with the remaining minimal phase gives 
the most compact profile compatible with the measured 
form factor \F\. In general, this is a good approximation, 



measurement errors may have larger impact (see, e.g., 
Refs. [301 H2])- Changing to wavelengths and removing 
the singularity at u) = uj' in Eq. ([6} result in |4TJH4"2] 
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1-(A'/A)2 

The inverse Fourier transform of Eq. ([s]) can be given by 

dX' 



Pn{z) 



|F(A')|cos(2^/A'-$(A')) 



A' 2 



(8) 



when considering only positive wavelengths (A > 0). The 
Eqs. Q and ([8| are the fundamental expressions for lon- 
gitudinal electron bunch profile reconstruction from spec- 
troscopic measurements and are used later in Sec. |III B| 



C. Statistical fluctuations of the radiation pulse 
energy in SASE FELs 

The amplification process in a SASE FEL starts from 
shot noise in the electron beam, passes the stage of ex- 
ponential amplification (s/L sat ~ 0.8), and finally enters 
the saturation regime (s/L sat ~ 1), where s/L sat is the 
normalized undulator length with L sat being the satura- 
tion length (see Fig. [3|. The FEL radiation pulse en- 
ergy exhibits shot-to-shot fluctuations, which are larger 
for shorter pulse durations. The maximum fluctuations 
are present at the end of the high-gain exponential gain 
regime. Radiation from a SASE FEL operating in the 
exponential gain regime possesses the properties of com- 
pletely chaotic polarized light [2Ql EI] • One consequence 
is that the probability distribution of the energy in the 
radiation pulse is given by the gamma distribution 
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where T{M) is the gamma function, M = 1 / '(a e I '(E)) 2 , 
and oe = \J {{E — (E)) 2 ) is the FEL radiation pulse en- 
ergy spread. The parameter M can be interpreted as the 
average number of "degrees of freedom" or "modes" in 
the radiation pulse (see, e.g., Refs. [2Ql[2T]). With knowl- 
edge of the coherence time r c , the FEL radiation pulse 
duration can be estimated as ~ M x r c . This estimate as- 
sumes a high degree of transverse coherence, which is true 
for the parameter space of well designed SASE FELs like 
FLASH. In the framework of a one-dimensional model, 
the maximum value of the coherence time 
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and the saturation length 
L ~ Au 
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are expressed in terms of the FEL parameter p and 
the number of the electrons N c — I / (epco) cooperating to 
the SASE FEL process [101 1211 S3 • Here, w/(2tt) = c/A 
is the frequency of the amplified wave, / is the electron 
bunch current, and A u is the undulator period. A suitable 
estimate for the parameter p comes from the observation 
that in the parameter range of SASE FELs operating 
in the extreme-ultraviolet and X-ray wavelength range, 
the number of field gain lengths to saturation is about 10 
(e.g., Ref. [2Q]). The FEL parameter p and the coherence 
time r c are related to the saturation length L sat as: 



K/L 



AL sat /(2^FcA u ) . (12) 



These simple physical considerations are confirmed with 
numerical simulations using the time-dependent simula- 
tion code FAST 05]. Here, we consider the model of a 
Gaussian longitudinal electron bunch profile, and trace 
FEL parameters for different values of the r.m.s. elec- 
tron bunch length at e . Figure 3(a) shows the evolution 



of the energy E in the radiation pulse, the fluctuations 
of the radiation pulse energy <je/(E), and the normal- 
ized r.m.s. FEL radiation pulse duration pwcr t p along 
the undulator. Qualitative observations are that the ra- 
diation pulse energy and pulse duration grow with in- 
creasing electron bunch length, and the maximum fluc- 
tuations of the radiation pulse energy and the minimum 
r.m.s. radiation pulse duration are observed at the end of 
the exponential gain regime (s/L sat < 0.8). The normal- 
ized values of these parameters, i.e., E/E sat , as/cr^ 3 *, 
and cr^p/crj 11 ™, exhibit nearly universal dependencies for 
a normalized electron bunch length pioo t ,c ^ 1 as it is 
shown in Fig. |3(b)[ This allows to derive a universal de- 
pendency between the r.m.s. electron bunch length a t ,e 
and the full width at half maximum (FWHM) radiation 
pulse duration T™ in = 2^2 ln2cr t ™ at the end of the 
exponential gain regime as a function of the number of 
modes M in the FEL radiation pulse (see Fig. Eh. For 
M > 2 and with reasonable practical accuracy we have: 
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(13) 



The minimum FWHM FEL radiation pulse duration ex- 



pressed in terms of the coherence time (cf. Eq. (12)) is 
T™ in ~ 0.7 x M X r c . Lengthening of the FEL radiation 
pulse occurs when the amplification process enters the 
saturation regime and happens due to two effects. The 
first effect is lasing toward saturation in the tails of the 
electron bunch, and the second effect is radiation pulse 
lengthening due to slippage effects, which is one radiation 
wavelength per undulator period. The effect of lasing in 
the tails gives the relative radiation pulse lengthening as 
it is illustrated in Fig. 3(a) At saturation (s/L S3i t ~ 1), 
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FIG. 3: Simulations on the evolution of FEL related parame- 
ters versus normalized undulator length s/L sat : (a) Evolution 
of the energy E in the radiation pulse (solid lines), the fluctua- 
tions of the radiation pulse energy <je/{E) (dashed lines), and 
the normalized r.m.s. FEL radiation pulse duration puuat tP 
(dotted lines) for normalized r.m.s. electron bunch lengths of 
puiot.c — 2 and 4 (blue and red), and (b) the normalized val- 
ues of E/E sat , o"E/o"E ax > an d o"t, P /o"™p n for the r.m.s. electron 
bunch lengths of pLda t , c = 2, 4, and 8 (blue, red, and green). 



Measurements of statistical properties of SASE FELs 
allow to estimate mean photon pulse durations and the 
lasing fraction of the FEL driving electron bunches, and 
this technique has been effectively used at FLASH and 
LCLS [U [U First, the FEL process is tuned to the 
maximum radiation pulse energy, which occurs in the 



saturation regime and beyond {s/L s 



> 



1), and then 



the pulse lengthening is a factor of about 1.4 with re- 
spect to the minimum pulse duration in t he e xponential 
gain regime (s/L sat < 0.8) given by Eq. (13 1, and it is 
increased up to a factor of 2 in the deep nonlinear regime 
(s/L sat > 1.2). It is also seen that the slippage effect is 
more pronounced for relative lengthening of short pulses. process is tuned to the end of the exponential regime 



the FEL radiation pulse energy is recorded at different 
positions along the undulators by applying orbit kicks 
in order to suppress lasing. The resulting "FEL gain 
curve" permits the estimation of the saturation length 
L sa t, and the FEL parameter p and the coherence time 
t c can be calculated by using Eq. (12 1. Then the FEL 
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FIG. 4: Normalized r.m.s. electron bunch length pcja t . a con- 
tributing to the FEL process (red dashed line) and the cor- 
responding normalized minimum FWHM radiation pulse du- 
ration pujT™ ln (blue solid line) at the end of the exponential 
gain regime (s/L sa t < 0.8) versus the number of modes M. 



at s/L s 



< 



0.8 where the FEL radiation power is re- 



duced by a factor of ~ 20 with respect to the saturation 
regime (see Fig. 3(a)). The fluctuations of the FEL ra- 



diation pulse energy reach their maximum value at this 
point, and the inverse squared value of radiation pulse 
energy fluctuations gives the number of the modes in the 
radiation pulse. The electron bunch length and the mini- 
mum photon pulse duration at the end of the exponential 
gain regime are then derived from Eq. (13). The relative 



photon pulse lengthening at saturation corresponds to a 
factor of about 1.4, and in the (deep) nonlinear regime, 
it depends on the electron bunch length (see Fig. |3(b)[ ) . 



III. EXPERIMENTAL RESULTS AT FLASH 

Comparative measurements of soft X-ray pulse dura- 
tions and electron bunch lengths were carried out for two 
accelerator and beam settings. In the first step, we mea- 
sured longitudinal electron bunch profiles determined by 
both time-domain longitudinal phase measurements and 
frequency-domain THz spectroscopy for electron bunch 
charges of 100 pC and 500 pC in order to verify the consis- 
tency of both methods. In the second step, we compared 
longitudinal electron bunch profiles determined by lon- 
gitudinal phase measurements with the soft X-ray pulse 
durations estimated by SASE FEL radiation pulse en- 
ergy statistics for electron bunch charges of 150 pC and 
500 pC. The relevant parameters are given in Table |TJ 



A. Longitudinal phase space measurements 

The longitudinal phase space measurements were per- 
formed using a 100 /jm thick scintillator (YAG:Ce) imag- 
ing screen, and the camera system was set up with a 



TABLE I: Parameters for the comparative measurements of 
FEL photon pulse durations and electron bunch lengths. The 
numbers in the brackets of the electron bunch charge and the 
FEL radiation saturation level are related to each other. 



Parameter 



Symbol Value Unit 



Electron beam energy 
FEL radiation wavelength 
Electron bunch charge 
Saturation level of FEL radiation 
Undulator period 



A 
Q 



660 MeV 

14.6 nm 

150 (500) pC 

30 (200) jitJ 

27.3 mm 



spatial resolution of better than 16 fim. The time cal- 
ibration was achieved by measuring the vertical elec- 
tron bunch centroid on the imaging screen during a 
TDS rf phase scan. By measuring the intrinsic vertical 
r.m.s. beam size, i.e., with the TDS rf power switched off, 
the r.m.s. time resolutions was estimated (see Sec 



IIA) 



for each measurement individually. For comparisons with 
soft X-ray pulse durations (see Sec. IV), 50 consecutive 



single-shot measurements of the longitudinal phase space 
were taken for each set of electron beam parameters with 
bunch charges of 150 pC and 500 pC. The longitudinal 
phase space and electron bunch current of typical single- 
shots, together with the single-shot r.m.s. bunch lengths, 
are shown in Fig. [5j The presented measurements were 
performed for one zero-crossing TDS rf phase, and cor- 
rections of a potential longitudinal-to- vertical correlation 
C y were neglected with S y 3> C y » 0. This is justified by 
the large shear parameter (see below) and the agreement 
in the comparative measurements between TDS and THz 
spectrometer presented in Fig. |6j where the latter is not 
sensitive to longitudinal-to-vertical correlations. 

The single-shot r.m.s. time resolution lZt,c reached un- 
precedented 8fs for the 150 pC case, which corresponds 
to a shear parameter of S y c~ l — 30 (significantly larger 
than any longitudinal-to-vertical correlation C v observed 
at FLASH) and an intrinsic vertical r.m.s. beam size of 
(T y o = 70 /im. For the 500 pC settings, we achieved an 
r.m.s. time resolutions of 13 fs, which is due to a larger 
intrinsic vertical beam size. The average r.m.s. electron 
bunch lengths, including statistical errors due to fluctu- 
ations of electron beam parameters and accelerator set- 
tings, of 50 consecutive single-shots are (E te ) = 41 ± 3 fs 
and 103±4 fs for the bunch charges of 150 pC and 500 pC, 
respectively. In order to take account of the unknown 
correlation C y , we admit a systematical uncertainty of 
±10 % in the r.m.s. electron bunch length measurements. 

The longitudinal phase space has a similar shape in 
both cases (Figs. 5(a) and |5(b)| ) with an overall negative 
energy chirp —dS/dt, i.e., electrons in the leading part of 
the bunch (on the left in Fig. [5| have less energy than 
those in the trailing part. However, both distributions 
show a distinct core region with a positive energy chirp, 
which is most probably generated by space charge forces. 
In general, the implementation of the new electron bunch 
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■£,«, = 41 ±2fs 



C_ - 1.5 




- 0.5 



-50 

t(fe) 

(a) Q = 150 pC. 




500 pC. 



FIG. 5: Single-shot longitudinal phase space and electron 
bunch current (dotted red lines) measurements using the 
TDS in combination with the energy spectrometer for bunch 
charges of (a) 150 pC and (b) 500 pC. The selected single- 
shots represent typical bunches out of 50 measurements, and 
the r.m.s. bunch length values (£t, e ) including measurement 
error in the legends correspond to the shown single-shots. 



compression scheme using the third third-harmonic rf lin- 
earizer system [251 |29J results in electron bunches with 
significantly more confined bunch current profiles with- 
out trailing tails of picosecond duration [21 [6] . 



B. Measurement of longitudinal form factors 

The principle of longitudinal electron beam diagnos- 
tics based on spectroscopy of coherent radiation has been 
described in Sec. |II B| The measurements of the longitu- 
dinal form factors were performed by spectroscopy of co- 
herent transition radiation generated from an aluminum 
coated silicon screen and using the THz spectrometer up- 
stream of the collimator (see Fig. [I]). Technical details 
about the THz spectrometer and the transport bcamlinc 
for THz radiation can be found in Refs. 38, 39, 46J. Sev- 
eral hundred single-shot THz spectra in the wavelength 
range 5-430 /xm were recorded for electron bunch charges 



of 100 pC and 500 pC. At the same time, single-shot lon- 
gitudinal phase space measurements using the TDS in 
combination with the energy spectrometer were carried 
out in order to compare and verify the reconstruction of 
the longitudinal electron bunch profiles from the spectro- 
scopic measurements. Diagnostics for FEL photon pulse 
durations were not available during these measurements. 

Figure [6] shows the measured mean moduli of the lon- 
gitudinal form factors, including \F\ derived from TDS 
measurements (red line), together with the reconstructed 
electron bunch currents in comparison with the cor- 
responding time-domain TDS measurements (insets of 
Fig. raj). The measurements of the longitudinal form fac- 
tors \F\ were done with two reflection grating sets ("THz: 



1 



THz: Short 
■ THz: Long 
TDS 



0.1 - 




20 40 60 



100 120 140 160 180 200 

A (yum) 

(a) Q = 100 pC. 



^ 0.1 




0.01 ^ 

~ 50 100 150 200 250 300 350 400 
A (yum) 
(b) Q = 500 pC. 

FIG. 6: Mean moduli of the longitudinal form factors mea- 
sured by THz spectroscopy of coherent transition radiation 
(green diamonds and blue squares) and derived from TDS 
measurements (solid red lines), generated by electron bunches 
with charges of (a) 100 pC and (b) 500 pC, respectively. The 
insets show single-shot electron bunch currents reconstructed 
from THz spectra (solid blue lines) in comparison with single- 
shot TDS measurements (dotted red lines), and the numbers 
in the brackets represent the calculated r.m.s. bunch lengths. 



Short" and "THz: Long") and averaged over ~ 300 sin- 
gle shots. The insets of Fig. [6] show single-shot elec- 
tron bunch currents reconstructed by means of an inverse 
Fourier transform using the Kramers-Kronig relations for 
the phase retrieval in comparison with direct measure- 
ments in the time-domain using the TDS. The integrals 
in Eqs. ([7]) and ^ for the electron bunch profile recon- 
struction were computed for single shots, representing 
the typical form factor out of the ~ 300 measurements, 
with a short-wavelength cut-off at 5 /Ltm and by applying 
an extrapolation with \F\ — > 1 for long wavelengths. 

The results of both experimental methods are in good 
agreement and show the same features within the elec- 
tron bunch currents. The deviations can be explained by 
uncertainties of both methods, including the unknown 
and omitted Blaschke phase in Eq. Q, and due to the 
fact that both diagnostics are not located at the same po- 
sition along the beamline. In fact, the collimator between 
the two positions (see Fig. [I]) generates longitudinal dis- 
persion, which changes the longitudinal electron bunch 
profile depending on the actual energy chirp of the elec- 
tron bunch and the settings of the collimator magnets. 
Nevertheless, the good overall agreement of the compara- 
tive measurement confirms the consistency of both high- 
resolution longitudinal electron beam diagnostics. 



C. Measurement of FEL pulse energy statistics 

The physical background of the FEL photon pulse du- 
ration estimations based on statistical properties of SASE 
FEL radiation has been described in Sec. Ill CI FLASH 
is equipped with the following set of detectors for single- 
shot measurements of the radiation pulse energy: gas 
monitor detectors, micro-channel plate (MCP) based de- 
tectors, photodiodes, and thermopiles [4TH5T)] . The de- 
tectors are installed at several positions along the pho- 
ton beamline (see also Fig. [I]). The MCP detector is 
installed in front of all other detectors and is used for 
precise measurements of the FEL photon pulse energy. 
The MCP measures the radiation scattered by a metal- 
lic mesh (Cu, Fe, and Au targets are being used) placed 
behind an aperture located downstream of the undulator 
magnets. The electronics of the MCP detector has an ex- 
cellent signal-to-noise ratio of about 100 and thus allows 
high-resolution measurements of the energy fluctuation. 

The measurement procedure was organized as follows. 
First, the SASE FEL process was tuned to the maxi- 
mum MCP signal at full undulator length of 27 m (six 
undulator modules). Then an orbit kick was applied by 
switching on steerer magnets after the fourth undulator 
module such that the FEL amplification process was sup- 
pressed in the last two undulator modules. The level of 
the FEL radiation pulse energy after four undulator mod- 
ules was about a factor of ~ 20 less compared to the level 
at saturation, and by taking into account the saturation 
length, which was estimated to be L sat = 22.5 ± 2.5 m for 
both electron bunch charges within the presented exper- 



imental conditions, the normalized undulator length was 
s/isat = 0.80 ±0.09 there. It has been shown in Sec. II C 



that this point corresponds to the end of the exponential 
gain regime. The aperture of the detector was adjusted 
such that all photons passed through the scattering tar- 
get of the detector. Thus, the MCP signal is proportional 
to the single-shot FEL radiation pulse energy. However, 
the fluctuations of electron beam and accelerator param- 
eters contribute to the fluctuations of the radiation pulse 
energies, while only fundamental SASE FEL fluctuations 
are essential. For this reason, a selection procedure had 
to be applied to the recorded data. The important elec- 
tron beam and accelerator parameters, e.g., bunch charge 
monitor readings, beam position monitor readings, read- 
ings of the bunch compression monitors, and read back 
values of rf parameters were recorded for each shot to- 
gether with the readings of the MCP detector. If electron 
beam or accelerator parameters deviated more than the 
prescribed threshold, the events were excluded from the 
data set. We note that the number of events in the data 
set after the selection procedure must be sufficiently large 
in order to provide an acceptable statistical accuracy. 

Figure [7] shows the probability distribution of the FEL 
radiation pulse energies for the electron bunch charges 
of 150 pC and 500 pC. The sets of raw data contained 
about 1400 (1700) single-shot measurements, and about 
800 (550) measurements remained after the data selec- 
tion procedure for 150 pC (500 pC). The selection was 
mainly achieved by removing large electron bunch charge 
fluctuations and rf phase drifts in the first accelerating 
structure (ACC1) upstream of the first bunch compres- 
sor. The data selection thresholds were chosen to remove 
clear correlations with the MCP readings. The measured 
number of modes amount to M — 2.8 ± 0.5 and 12 ± 2 
for the settings with 150 pC and 500 pC, respectively, and 
the quoted errors are based on uncertainties of the data 
selection procedure. If a Gaussian shape for the lasing 
fraction of the electron bunch is assumed, we can apply 



Eq. ( 13 1 in order to determine the r.m.s. electron bunch 



length that contributes to lasing and get <Jt,c — 23 ± 4fs 
and 96 ± 18 fs for the electron bunch charge of 150 pC and 
500 pC, respectively. The minimum FWHM photon pulse 
durations T™ 111 at the end of the exponential gain regime 
are estimated to be these obtained r.m.s. electron bunch 
lengths (see Eq. (pj)), i.e., T' nin = 23±4fs and 96±18fs, 



respectively. The quoted errors are purely statistical and 
assume a Gaussian shape for the lasing fraction of the 
electron bunch. Significant systematical errors can ap- 
pear due to large deviations from the assumed longitu- 
dinal electron bunch profile and time-dependent electron 
bunch properties such as the transverse slice emittance. 

In parallel with statistical measurements, FEL photon 
spectra were recorded using the soft X-ray (SXR) spec- 
trometer in the photon beamline. We note that the spec- 
trum is simply a Fourier transform of the temporal struc- 
ture, and the average number of spikes (modes) in the 
time-domain should be about the same as the number of 
spikes in the spectral-domain. The single-shot soft X-ray 
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0.5 1 1.5 2 2.5 



E/(E) 
(b) Q = 500 pC. 

FIG. 7: The probability distributions p(E) of the FEL radia- 
tion pulse energy E of soft X-ray pulses at the end of the ex- 
ponential gain regime for electron bunch charges of (a) 150 pC 
and (b) 500 pC. The histograms (blue bars) represent the data 
with the relative radiation pulse energy spread oe/(E), and 
the solid curve (red lines) show the gamma distribution using 
M — 2.8 and 12 (modes) for 150 pC and 500 pC, respectively. 



spectra were obtained by using a plane grating monochro- 
mator, operating in a spectrographic mode with a resolv- 
ing power A/ A A > 10000 [ST]. We observed qualitative 
agreement between the spectral measurements (number 
of spikes) and the number of modes measured by FEL 
pulse energy statistics (see Fig. [8]). The detailed analy- 
sis of the spectral measurements, taking into account the 
energy chirp of the electron bunches, is beyond the scope 
of this paper and will be reported elsewhere. 

Based on the theore tical m odel described in Sec. Ill CI 
and according to Fig. |3(b)| the expected FEL photon 
pulse durations at saturation and beyond can be extrap- 
olated by the minimum FWHM photon pulse durations 
measured at the end of the exponential gain regime. In 
the nonlinear regime sit s I L sa i ~ 1.2 (sec Fig. [3(b)] ) and 
by using the measured minimum FWHM photon pulse 




A (urn) 
(b) Q = 500 pC. 

FIG. 8: Single-shot (red, green, and blue solid lines) and mean 
(dashed black lines) spectra of soft X-ray pulses measured at 
the end of the exponential gain regime using the SXR spec- 
trometer for the settings with (a) 150 pC and (b) 500 pC. The 
mean spectra represent an average over ~ 2500 single shots. 

durations T™, we would expect FWHM FEL photon 
pulse durations of T p = 43 ± 8fs (lengthening of ~ 1.9) 
and T p = 169±32 fs (lengthening of ~ 1.75) for the beam 
settings with 150 pC and 500 pC, respectively. For these 
expected FEL photon pulse durations and by considering 
the measured average FEL pulse energies (see Table |T] of 
30 /iJ (150 pC) and 200 [iJ (500 pC), the peak powers re- 
sult in 0.7 ± 0.1 GW and 1.1 ± 0.2 GW, respectively. 



IV. SOFT X-RAY PULSE DURATIONS VERSUS 
ELECTRON BUNCH LENGTHS 

The average r.m.s. electron bunch lengths determined 
by time-domain longitudinal phase space diagnostics are 
(S t o ) = 41 ± 3fs and 103 ± 4fs for the bunch charges 
of 150 pC and 500 pC, respectively. The corresponding 
minimum FWHM soft X-ray pulse durations at the end 



10 



1.5 



H 



0.5 





- - - TDS measurement - 
. — Electron bunch: cr^ c 

V T~) 1_ — 1, — ,„ — mill 

* rnoton pulse: a t p 




F -V 
1" "\".* 


* • / ■ f 

* */•" T 
i »' / ; -1." 

* • / « '7* 

* * / .* 

* * / • ■1." 

* •* / * 

* * # * ■'/»" 


t : V. *» 

V ■ \ * * 

* A- *» • 

■\*. * \ *• * 
"\" * \ *. * 
•V' *• \ ** * 
*v« * \** * 
•V ** \.** * 







-100 



-50 





t(fe) 



50 



100 



FIG. 9: Single-shot longitudinal electron bunch profile mea- 
surement using the TDS for a bunch charge of 150 pC (dashed 
red line). The selected single-shot represents a typical bunch 
out of 50 measurements (cf. Fig. 5(a) I with an r.m.s. bunch 
length of Et, c = 41 ± 2fs. The blue solid (dotted) line shows 
a Gaussian with a — 23 ± 4fs representing the width a t , c of 
a potential lasing part of the electron bunch. The green solid 
(dotted) line shows a Gaussian with a = (23±4)/(2V21n2) fs 
representing the corresponding soft X-ray (photon) pulse du- 
ration (T™ 1 . The dotted lines describe the error bands (±). 



of the exponential gain regime estimated by SASE FEL 
pulse energy statistics are T" in = 23±4fs and 96±18fs, 
respectively for the same electron bunch charges. 

As described in Sec. |II C| for Gaussian longitudinal 
profiles, the r.m.s. electron bunch length cr tiC that con- 
tributes to lasing is related to the minimum FWHM 



FEL photon pulse duration T t 
where T, min 



mm 
P 



by (T t ,, 



o^min 
P 



(see 



Eq. (13 1), where T™ in = 2 V2 In 2 a™ n with the corre- 
sponding minimum r.m.s. FEL photon pulse duration 
a - ™™. The comparison of £t. c and T™" 1 (at.c) shows that 
the r.m.s. electron bunch lengths measured by longitu- 
dinal phase space diagnostics are larger than estimated 
by SASE FEL pulse energy statistics, which is more sig- 
nificant for the bunch charge of 150 pC. Figure [9] shows 
a measured single-shot longitudinal electron bunch pro- 
file (I in kA) for 150 pC and two Gaussians (I in arb. 
units) with standard deviations a — <Jt,e — T^ un , which 
represents a potential lasing part of the electron bunch, 



and 



>t, v 



= T™ n /(2v / 21n2), representing the cor- 



responding soft X-ray pulse. The amplitudes and peak 
positions of the Gaussians were chosen to fit the core 
region of the measured longitudinal electron bunch pro- 
file. The deviations of the estimated minimum FWHM 
soft X-ray pulse duration T™ m , i.e., the lasing fraction 
of the electron bunch tr t e , from the measured r.m.s. elec- 
tron bunch length E tje can be explained most likely by 
the non-Gaussian longitudinal electron beam profile with 
low currents in the tails that do not contribute to lasing 
but to the r.m.s. electron bunch length value. 

In cases of non-Gaussian longitudinal electron bunch 
profiles, e.g., an electron bunch compression mode with 



a short leading spike (e.g., Refs. [21 E]), a Gaussian-like 
core region needs to be identified in order to determine 
the pulse durations. For instance, the FWHM bunch 
length value, which is not that sensitive to tails in the 
electron bunch, can be determined and used to calculate 
an effective r.m.s. bunch length for comparison with the 
minimum FWHM FEL photon pulse duration. For the 
settings with 150 pC, the FWHM electron bunch length is 
70±17fs which corresponds to an effective r.m.s. electron 
bunch length of 30 ± 7 fs (cf. with T™ in = 23 ± 4 fs). 

For a more precise identification of the lasing part 
within the electron bunches, simultaneous information 
of other parameters, e.g., the slice emittance, or time- 
domain measurements of the FEL photon pulse profiles 
(e.g., Refs. [THlUnilSS]) are required. However, the longi- 
tudinal electron beam diagnostics discussed here provide 
robust estimates on the upper limits of the expected FEL 
photon pulse durations and allows their monitoring. 



V. SUMMARY AND CONCLUSIONS 

Information on the FEL photon pulse durations is im- 
portant for time-resolved and high intensity experiments 
in photon science. However, direct measurement of fem- 
tosecond FEL photon pulse durations in the extreme- 
ultraviolet and X-ray wavelength range is a tremendous 
challenge, and well established methods are not yet avail- 
able. The theoretical considerations, simulations, and ex- 
perimental data presented in this paper show that high- 
resolution electron beam diagnostics can be utilized to 
provide reasonable constraints on expected FEL photon 
pulse durations. We discussed two different longitudinal 
electron beam diagnostics with single-shot capability and 
femtosecond accuracy, and presented consistent single- 
shot measurements of longitudinal electron bunch pro- 
files with good agreement. The electron beam measure- 
ments were performed in the frequency-domain by THz 
spectroscopy of coherent transition radiation and in the 
time-domain by longitudinal phase space measurements 
using a TDS in combination with a magnetic energy spec- 
trometer. In the latter case, an unprecedented single- 
shot r.m.s. time resolution of 8 fs has been achieved. The 
results from the longitudinal phase space measurements 
were compared with soft X-ray pulse durations estimated 
by SASE FEL pulse energy statistics, and the consistent 
results show that r.m.s. electron bunch lengths, deter- 
mined by longitudinal electron beam diagnostics, provide 
upper limits on the expected FWHM FEL photon pulse 
durations. Theoretical considerations and FEL simula- 
tions verify this experimental observation and facilitate 
the extrapolation of the expected FEL photon pulse du- 
rations in the saturation regime and beyond based on the 
measurements in exponential gain regime. In addition, 
we demonstrated the generation of soft X-ray pulses with 
a peak radiation power in the GW level and with dura- 
tions well below 50 fs (FWHM) after implementation of 
the new uniform electron bunch compression scheme us- 
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